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REFLECTIVITY OPTIMIZATION FOR MULTILAYER STACKS 
FIELD OF THE INVENTION 

[001] The present invention relates to lithography optimization techniques, and 
more particularly to a method for obtaining an optimal reflectivity value for 
complex multilayer stacks. 

BACKGROUND OF THE INVENTION 

[002] Photolithography is a lithographic technique used to transfer the design of 
circuit paths and electronic elements of a chip onto a wafer's surface. A 
photomask is created with the design for each layer of the board or wafer (chip). 
The board or wafer is coated with a light-sensitive film (photoresist) that is 
hardened when exposed to light shining through the photomask. The board or 
wafer is then exposed to an acid bath (wet processing) or hot ions (dry 
processing), and the unhardened areas are etched away. 

[003] Ideally, the photoresist pattern produced by the photolithography process 
and the substrate pattern produced by the subsequent etch process would 
precisely duplicate the pattern on the photomask. For a variety of reasons, 
however, the photoresist pattern remaining after the resist develop step may vary 
from the pattern of the photomask significantly. 

[004] For example, many of the integrated circuit elements formed on the 
wafer*s surface comprise multiple layers of thin film, interferences from which can 
result in critical dimension (CD) variation. The interferences are mainly caused 
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by the reflectivity from resist top and bottom interfaces, and can lead to the 
existence of standing waves and CD swing curves. 

[005] FIG. 1 is a cross-sectional view of an incident of a light beam shown on a 
film stack shown during lithography. The light 10 shines through the resist 12, 
which is also referred to as the ambient, to the film stack 14, A portion of the 
light 18 passes through the film stack 14 into the substrate 16, another portion of 
the light 20 is reflected from the surface of the stack 14, and yet another portion 
of light 22 is reflected at boundaries between the layers in the stack 14. Optical 
functions, such as reflectivity of layer boundaries, can be computed by a set of 
parameters that include light wavelength, thickness of the stack layers, and a 
complex index of refractivity for the ambient, substrate, and each film layer. 

[006] The optimization for the lithography application typically involves the 
minimization or maximization of reflectivity. The optimization can be carried out 
by finding values for the parameters that cause a cost function of reflectivity to 
yield an optimal value. It is difficult or very time consuming to obtain an optimal 
value for complex multilayer stacks 14. 

[007] Software programs are available that simulate lithography process steps 
and parameters. PROLITH™ by KLA-Tencor Corporation of San Jose, California 
is an example of a standard lithography simulation program. PROLITH is capable 
of aerial image and three-dimensional resist image predictions. Although current 
lithography simulation programs, such as PROLITH, help lithographers reduce 
process development and process optimization times, such programs have 
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limitations. PROLITH, for example, performs the reflectivity minimization 
calculation by fixing two variables and finding the minimum of the function using 
the remaining variables. Two other variables are fixed, and the process is 
repeated for each layer. Thus, using conventional lithography simulation 
programs is a very iterative and time consuming process. In addition, results are 
not usually very accurate. 

[008] Accordingly, what is needed is an improved method for obtaining an 
optimal reflectivity value for complex multilayer stacks. The present invention 
addresses such a need. 

SUMMARY OF THE INVENTION 

[009] The present invention provides a method for obtaining an optimal 
reflectivity value for complex multilayer stacks. Aspects of the present invention 
include generating a model of a multilayer stack and parameterizing each layer 
by a thickness and an index of refraction; allowing a user to input values for the 
parameters; calculating an extrema for a cost function of reflectivity R using the 
input parameter values; calculating sensitivity values S for the extrema points; 
and obtaining an optimal value by calculating a cost function R + S. 

[010] According to the method and system disclosed herein, the present 
invention provides a method for obtaining an optimal reflectivity value for 
complex multilayer stacks that is capable of perfomning the calculation using 
multiple variables. In addition, the method is faster and more accurate than prior 
methods. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[011] FIG. 1 is a cross-sectional view of an incident of a light beam shown on a 
film stack shown during lithography. 

[012] FIG. 2 is a flow diagram illustrating the process of obtaining an optimal 
reflectivity value for a complex multilayer stack in accordance with a preferred 
embodiment of the present invention. 

[013] FIG. 3 is a diagram illustrating an example modeled stack. 

[014] FIG. 4 is a diagram illustrating the parameters of the multilayer stack. 

[015] FIG. 5 is a diagram showing a parameter value input screen for the 
reflectivity optimization. 

[016] FIG. 6 is a graph illustrating example sensitivity values for a cost function. 

[017] FIGS. 7A and 7B are diagrams showing an example reflectivity simulation 
and experimental CD swing curve as a function of resist thickness. 

DETAILED DESCRIPTION OF THE INVENTION 

[018] The present invention relates to obtaining an optimal reflectivity value for 
complex multilayer stacks. The following description is presented to enable one 
of ordinary skill in the art to make and use the invention and is provided in the 
context of a patent application and its requirements. Various modifications to the 
preferred embodiments and the generic principles and features described herein 
will be readily apparent to those skilled in the art. Thus, the present invention is 
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not intended to be limited to the embodiments shown, but is to be accorded the 
widest scope consistent with the principles and features described herein. 

[019] The present invention provides a computer-implemented method for 
obtaining an optimal reflectivity value for a complex multilayer stack during 
lithography simulation. In a preferred embodiment, the present invention is 
implemented as a program that is separate, or part, of a lithography simulation 
program, and is executed on a conventional workstation, PC. or server. 

[020] FIG. 2 is a flow diagram illustrating the process of obtaining an optimal 
reflectivity value for a complex multilayer stack in accordance with a preferred 
embodiment of the present invention. The process begins in step 100 by 
generating a model of a multilayer stack by adding layers of material on top of a 
substrate material, and parameterizing each layer by a thickness and an index of 
refraction in order to compute the optical functions. 

[021] FIG. 3 is a diagram illustrating an example modeled stack 50. A stack 
may comprise N layers. However, in the example shown, stack 50 comprises a 
total of 6 layers: a top ambient resist layer 0 followed by four layers of materials 
(1, 2, 3, and 4) patterned over a substrate layer 5. For purposes of example, the 
materials comprising layers 1-4 are Si02, SiOn, ILD, and Si3N4, respectively. All 
layers, ambient and substrate are supposed to be homogeneous, with abrupt and 
parallel border interfaces. Well-known Fresnel equations are calculated at 
normally incident plane wave with monochromatic illumination. 
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[022] FIG. 4 is a diagram illustrating the parameters of the multilayer stack 50. 
Assume that the multilayer stack comprises of N layers (L = t, 2, 3, N). Each 
layer 52 is parameterized by a thickness 54 and an index of refraction 56 (which 
includes a real and an imaginary number), including the ambient resist layer 0 
and the substrate layer N+1. The/^ layer has thickness ofy, and complex index of 
refraction nj = nj - i kj, as shown. The ambient and substrate, which are above 
and below the multilayer stack, have complex indexes of refraction: no = no - i ko 
and Hn+i = nN+i - i kN+i, respectively. 

[023] The reflectivity 55 at an interface between two layers 52 can be defined as 
a cost function. For example, the reflectivity 55 at the interface (between the (y 
- and f layers) is Ry, is a function of 3{N - y + 1) + 4 parameters, which are: 
Hy-t, nj ... n/v, H/v+t; /cy.t, kj ... /c/v, /c/v+t; dy, ofy+y ... of/v. For lithography applications, 
the resist is usually on the top layer of multilayer stacks so that the cost functions 
(Ri and R2) are the reflectivity from the top and bottom of the resist. They can 
have 3A/ + 4 and 3 A/ + 1 independent parameters, respectively. 

[024] Based on practical problems, the cost functions can be evaluated and 
calculated by fixing some parameters and varying others. Local extrema (minima 
or maxima) of the cost functions of reflectivity for these multidimensional cases 
with defined feasible ranges may be solved by a well-known variable 
neighborhood method. 

[025] Referring again to FIG. 2, after the stack 50 has been modeled, in step 
102, the user inputs values for the parameters, including range and step values. 
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FIG. 5 is a diagram showing a parameter value input screen for the reflectivity 
optimization. The user may first enter values for the thickness 54 and the 
complex indexes of refraction (n and k) 56 and 58 for each layer 52. The values 
include a current starting point 60, a minimum values 62, and a maximum value 
64. The user also chooses which of the parameters will be independent 
variables with feasible intervals by clicking the "Vary" button 59 and inputting a 
step value 66. Those parameters that are not designated as varying will be fixed. 
According to one aspect of the present invention, the user may vary any number 
of parameters. The user also chooses a particular layer for the reflectivity 
calculation (not shown). 

[026] Referring again to FIG. 2, after the user has input parameter values, in 
step 104, the extrema (min and max) for the cost function of reflectivity R is 
calculated for the top of the selected layer Interface using the input parameter 
values. 

[027] Since the variability in film stacks in semiconductor processing can lead to 
a variation of the parameters, the sensitivity of reflectivity to the parameters has 
to be taken into consideration in order to find an optimal condition. Therefore, in 
step 106, all sensitivity values S for the extrema points are calculated. The 
sensitivity is defined as S = (Max R - Min R) with all varied parameters, 
which represents how reflectivity values react to the variation of variables. The 
larger sensitivity, the higher magnitude of change of reflectivity should be. 
Smaller sensitivity values are preferred since the variation of thickness and 
refractive coefficients of materials results In smaller fluctuation of reflectivity. 
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[028] FIG. 6 is a graph illustrating example sensitivity values for a cost function. 
The global minimum value for the function is point A, but the responding 
sensitivity is high. The global minimum value for the reflectivity is point A, but the 
responding sensitivity in that point is high. Point B is an optimal point because 
the sensitivity value is much smaller than the value at point A though 
the reflectivity value is slight higher. In order to obtain the most favorable 
solution, another cost function, R + a S, is computed in a preferred embodiment, 
in which a is a weighted parameter for sensitivity in the cost 
function. The larger alpha value, the more important the sensitivity is. 

[029] Referring again to FIG. 2, after the sensitivity values are calculated, in 
step 108, a value for the cost function R + a S is calculated, from which an 
optimal value is obtained in step 110. The choice of the best optimal point also 
may depend upon other process considerations such as etch requirements. 

[030] The simulation of reflectivity for a multilayer film stack was validated using 
experimental CD swing results since the resist CD swing curve is controlled by 
the refection of light at the resist boundaries. For the experiment, a film stack 
comprising the following layers was used: 
Ambient: air 

Layer # 1 : DP630, thickness = 290 - 600 nm 
Layer # 2: SiON, thickness = 120 nm 
Layer # 3: Black Diamond, thickness = 770 nm 
Layer # 4: SiCO, thickness = 50 nm 
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Substrate: Cu 



[031] In addition, the illumination wavelength Is 248 nm. The reflectivity 
simulation and experimental CD curve as a function of resist thickness are shown 
in FIGS. 7A and 7B, respectively. The Standing Wave Ratio (SWR) is defined 
as (CD maximum - CD minimum)/(CD maximum + CD minimum) * 100. A lower 
SWR means lower developed resist CD variation due to resist thickness and 
substrate variation. The agreement between reflectivity simulation and 
experiment is excellent. The simulated reflectivity minima correspond to the 
bright field experimental developed resist CD maximums and the simulated 
reflectivity maxims correspond to the experimental developed resist CD minima. 
The extema locations are quantified in Table 1 . 

[032] Table 1 . Resist values at extremum points 



Experiment (nm) 


Simulation (nm) 


342 


342 


378 


378 


413 


413 


448 


448 


484 


485 


519 


518 
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[033] A method for obtaining an optimal reflectivity value for complex multilayer 
stacks has been disclosed that is capable of performing the reflectivity calculation 
using multiple variables at once. In addition, the method is faster and more 
accurate than prior methods. 

[034] The present invention has been described in accordance with the 
embodiments shown, and one of ordinary skill in the art will readily recognize that 
there could be variations to the embodiments, and any variations would be within 
the spirit and scope of the present invention. Accordingly, many modifications 
may be made by one of ordinary skill in the art without departing from the spirit 
and scope of the appended claims. 
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